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Abstract. While the inspiral and ring-down stages of the binary black-hole coalescence 
are well-modelled by analytical approximation methods in general relativity, the recent 
progress in numerical relativity has enabled us to compute accurate waveforms from the 
merger stage also. This has an important impact on the search for gravitational waves 
from binary black holes. 'Complete' binary black-hole waveforms can now be produced 
by matching post-Newtonian waveforms with those computed by numerical relativity, 
which can be parametrised to produce analytical waveform templates. The 'complete' 
waveforms can also be used to estimate the efficiency of different search methods aiming 
to detect signals from black-hole coalescences. This paper summarises some recent efforts 
in this direction. 



1. Introduction 

Coalescing black-hole binaries are among the most promising sources of gravitational 
waves for ground-based detectors like LIGO and Virgo, and the planned space-borne 
detector LISA. The evolution of binary black holes is conventionally split into three stages: 
inspiral, merger and ring down. In the inspiral stage, the two compact objects, driven 
by radiation reaction, move in quasi-circular orbits. Eventually approaching the ultra- 
relativistic regime, the two bodies merge to form a single excited Kerr black hole. In 
the ring-down stage, the excited black hole loses its energy by gravitational-wave emission 
and settles into a Kerr black hole. Gravitational waveforms from the inspiral and ring- 
down stages can be accurately computed by approximation/perturbation techniques in 
general relativity [XJ [2j [3] . The recent progress in numerical relativity [4j [5j E] has enabled 
us to model also the non-perturbative merger phase of the coalescence of binary black 

holes ia e ei EBini iia nana nana nzi- 

While the current gravitational-wave searches look for each stage of the binary black- 
hole coalescence separately (see, for example, [H1IH]), combining the results from analytical 
and numerical relativity enables us to coherently search for all the three stages using a 
single template family. This coherent search is significantly more sensitive than the current 
searches over certain mass ranges (see Section [1]). This search has added advantages: 
including all the three stages adds more 'structure' to the template waveform, resulting 
in a potential reduction of false alarms. The additional structure and the improved signal- 
to-noise ratio also results in an improved estimation of the parameters of the binary, which 
is particularly important for LISA data analysis. As LISA data will contain a 'cocktail' 
of many strong binary signals, these will have to be subtracted from the data in order 
to analyse other signals. Improved parameter estimation can also have a tremendous 
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impact in cosmology. Since many of the supermassive black-hole mergers are likely to 
have electromagnetic counterparts, it is possible to constrain the values of cosmological 
parameters by combining the gravitational- wave and electromagnetic observations |20j . 
In particular, using the distance-redshift relation from many binary black-hole 'standard 
sirens\ LISA might be able to put interesting constraints on the equation of state of the 
dark energy The error bars on this depend on how accurately the 'red-shifted' mass of 
the source and the luminosity distance are estimated, and how well the host galaxy of the 
electromagnetic counterpart is identified. The improved parameter estimation might help 
to tighten these constraints. 

Several authors have proposed different ways of computing gravitational- wave templates 
containing all the three stages of the binary black-hole coalescence [331 [2U [3S1 [25] . In 
particular, Ref. [25] proposed a phenomenological parametrisation for non-spinning binary 
black-hole waveforms. These waveforms are explicit functions of the physical parameters of 
the system and exhibited very high overlaps with the 'target signals'. The target signals were 
constructed by matching numerical-relativity waveforms with post-Newtonian waveforms in 
appropriate matching regions. These waveforms contain all the three stages of the binary 
black-hole coalescence. Hence, they can also be used to estimate the efficiency of different 
(template-based and other) search methods used to detect signals from binary black-hole 
coalescences. 

This paper provides an overview of the application of the results from analytical and 
numerical calculations of binary black-hole waveforms into gravitational-wave data analysis. 
Section [2] describes how 'complete' binary black-hole coalescence waveforms (so-called hybrid 
waveforms) from non-spinning binaries can be constructed by matching post-Newtonian and 
numerical- relativity waveforms. Here we also investigate the robustness of the matching 
procedure by studying the mismatch between hybrid waveforms constructed using different 
matching regions. Section [3] introduces an analytical two-parameter family of non-spinning 
waveforms having very good overlaps with the hybrid waveforms. Section [4] shows how the 
hybrid waveforms can be used to estimate the sensitivity of different searches. Here we 
compare the fitting factor and faithfulness of different template-based searches using the 
hybrid waveforms as target signals. Finally, Section [5] provides a summary and future plans. 



2. Constructing 'complete' binary black-hole coalescence signals 

Although numerical relativity (NR) is able to compute gravitational waveforms containing 
all the three stages of the binary black-hole coalescence, the numerical simulations are 
heavily limited by the computational resources. But, the post-Newtonian (PN) formalism is 
known to work very well in the early inspiral. Thus, complete binary black-hole coalescence 
waveforms (hybrid waveforms) can be constructed by matching PN and NR waveforms in 
an appropriate matching region. Different authors have studied the consistency of the non- 
spinning PN waveforms with NR waveforms. See Refs. [27l Ell [Ml Ell [26j [30j El [25l [26] for 
some of the recent work. 

The time-domain waveform h+ >x (t,fji) from a particular system is parametrised by a 
set of 'extrinsic parameters' fi — {v?o,io} 5 where <po is the initial phase and to is the start 
time of the waveform. The PN and NR waveforms are matched by minimising the integrated 
squared difference, 5, between them in the matching region t\ <t<t%: 

dt. (1) 

i=+,X " L1 



K {t,n)-ah l (t,fj,) 
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Figure 1. Top-left panel shows an example hybrid waveform (green) constructed 
by matching an equal-mass NR waveform (red) computed by the Jena group with a 
3.5PN restricted PN waveform (black dashed). The two black vertical lines indicate the 
matching region (— 280A/ < t < —206M) employed. The rest of the panels show the 
mismatch between hybrid waveforms h n and ?t n _i, where a subscript n means that the 
hybrid waveform is constructed by matching the NR and PN waveforms at the n th cycle 
of the NR, waveform. The horizontal axis reports the start time, ti, of the nth cycle (in 
units of M). The mismatch is computed using Initial LIGO (top-right), Advanced LIGO 
(bottom-right) and Virgo (bottom-left) noise spectra. Total mass of the waveforms (in 
units of Mq) is shown in the legends. A mismatch of 3% is marked with a dashed 
horizontal line. 



The minimisation is carried out over the extrinsic parameters \x of the PN waveform and an 
amplitude scaling factor a [J] while keeping the 'intrinsic parameters' (the two component 
masses) of both the PN and NR waveforms the same. The hybrid waveforms are then 
produced by combining the 'best-matched' PN waveforms with the NR waveforms in the 
following way: 

(*, n) = a r(t) h N + R x (t, /*) + (1 - r(t)) h™ x (t, Mo) (2) 
where fiQ and oq denote the values of fi and a for which S is minimised, and r is a weighting 

X The NR waveforms used for this work contains only the I = 2, m = ±2 modes, and the PN waveforms 
are computed in the restricted PN approximation. Since the PN corrections to the amplitude of the PN 
waveforms are ignored, this will introduce an error of ~ 8% in the amplitude of the hybrid waveforms, 
and hence in the horizon distances reported in Figure [3] but not in the calculation of fitting factor and 
faithfulness. 
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function, denned as 



if t < t t 



r(t) = { if * 1 <*<t 2 (3) 

1 if < 2 < t. 

It is expected that the early inspiral is better-modelled by PN waveforms, as the early- 
inspiral part of the NR waveforms are prone to larger errors. But the PN waveforms become 
less accurate at the late inspiral. Thus, it is important to choose a matching region where 
both waveforms are accurate. This also means that the hybrid waveforms constructed using 
two very different matching regions can potentially be quite different. If, on the other hand, 
we are able to show that these differences are not very significant for data-analysis purpose, 
this is an indication that our analysis is not heavily dependent on the choice of the matching 
region. 

As a test of the robustness of the matching procedure, we compute the mismatch 
between two hybrid waveforms h n and h n -i, defined as 



MM(h n , h n -i) = 1 - max to 



4 Re 



M/)^-i(/)e i2w/t °d/ 



(4) 



s h {f) 

where Sh{f) is the one-sided power spectral density of the detector noise. The subscript n 
on the hybrid waveform h means that the hybrid waveform is constructed by matching PN 
and NR waveforms at the nth cycle of the NR waveform. 

The top-left panel of Figure [l] shows an example set of the PN, NR and the hybrid 
waveforms. The hybrid waveform is constructed by matching an equal-mass NR waveform 
computed by the Jena group, reported in Ref. [30j , with a restricted 3.5PN TaylorTl |32j 
waveform. Other panels in Figure [T] show the mismatch between the hybrid waveforms h n 
and h n -i, computed using three different noise spectra. If we take 3% as the maximum 
allowed mismatch between hybrid waveforms h n and h n —i, this preliminary exercise suggests 
that any matching region before t\ = — 150Af is robust for constructing hybrid waveforms 
using the equal-mass NR waveforms considered here. This will be studied in detail in a 
forthcoming work [33 . 



3. Templates for binary black-hole coalescence 

The hybrid waveforms constructed in the previous section can be parametrised in terms of 
the two physical parameters of the binary, thus producing analytical waveform templates. 
These analytical waveforms can be used to construct template banks for matched-filter 
searches; thus avoiding the computational cost of generating hybrid waveforms at each grid 
point in the parameter space. Ref. |25j proposed a family of Fourier domain templates of 
the form: 

u(f)=A eS (f)e^f\ (5) 
where the effective amplitude and phase are expressed as: 

f (///mcrg)" 7/6 if /< /mcrg 
Aes(f) = C < (/// mcrg r 2/3 if / mcrg < / < / ring 
I wC(f, /ring, a) if / r ing < / < /cut, 
1 7 

tfeff(/) = 27Tft + <p + - V(zfc V 2 + VkV + Zk) (7rM/)( fc - 5 )/ 3 . (6) 

77 * — ' 

' fe=0 
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In the above expressions, C is a numerical constant whose value depends on the sky-location 

and orientation of the binary as well as its physical parameters. For optimally-located and- 

M s/6 f-T/e re- 
oriented binaries, C = — d7r 2/3° rs y ^4 • to is the time of arrival of the signal at the detector, 

ifct the initial phase, £(/, /ring, o") = (5^) (/_/ °)2+o- 2 /4 a Lorentzian function with width <r 
centered around the frequency / r i ng , w a normalisation constant chosen so as to make A c g(f) 
continuous across the 'transition' frequency /ring, and / mer g is the frequency at which the 
power-law changes from /~ 7 / 6 to /~ 2 / 3 . The phenomenological parameters / mcrg , / r i ng , a 
and / cut are written in terms of the total mass M and symmetric mass ratio 77 of the binary 
as 

7rM/ mcrg = a i] 2 + b 77 + cq , 
tM / ring = ai ?7 2 + 61 77 + ci , 
7rMer = a 2 rj 2 + b 2 t] + c 2 , 

nMf cut = a 3 ?7 2 + 6377 + C3. (7) 

The coefficients aj,bj,Cj, j — 0...3 and Xk,yk,Zk, k — 0,2,3,4,6,7 are unique for a 
given family of hybrid waveforms. The coefficients corresponding to the hybrid waveforms 
considered here are tabulated in Table [T] These are computed from 7 hybrid waveforms in 
the range 0.25 > r/ > 0.16 produced by matching the numerical waveforms (referred to as 
the 'long NR waveforms' in Ref. [3S]) produced by the Jena group with restricted 3.5PN 
TaylorTl waveforms. 
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1.8046 xlO 2 


1p7 












Table 1. Coefficients describing the amplitude and phase of the phenomenological 
waveforms. See Eqs. (16b and JTB . 



4. Assessing the efficiency of different searches 

The hybrid waveforms constructed in Section [2] can be used to estimate the efficiency of 
different searches (including the one proposed in the previous section) in detecting signals 
from binary black-hole coalescences [|] Ideally, this should be done by injecting a large 
number of hybrid waveforms from different binaries into the detector data and by estimating 

§ The obvious assumption involved is that the hybrid waveforms are sufficiently close to the signals produced 
by nature. This is greatly dependent on the systematic errors in the hybrid waveforms. 



G 




Figure 2. Fitting factors (left plots) and faithfulness (right plots) of the PN inspiral, ring 
down and coalescence templates. Dots correspond to PN inspiral templates, triangles to 
ring-down templates, and squares to coalescence templates. The overlaps are computed 
using three different noise spectra — Initial LIGO (thin lines), Virgo (thicker lines) and 
Advanced LIGO (thickest lines). Horizontal axes report the total mass of the binary. 



the fraction of the injections detected by each search. Different search groups have already 
started work in this direction [M]. This section presents a simple strategy to estimate 
the efficiency of different template families in detecting signals from binary black-hole 
coalescences. The template families being considered here are (i) restricted 3.5PN TaylorTl 
inspiral templates truncated at the maximum binding-energy circular orbit |35j (ii) black- 
hole ring-down templates proposed in Ref. [35], and (iii) black-hole coalescence templates 
described in the previous section. 

We compute the fitting factors [37] and faithfulness [38] of different template families 
with the hybrid waveforms. Fitting factor is the overlap of a template waveform with the 
target signal maximised over both the intrinsic (M and 77) and the extrinsic (to and ipo) 
parameters of the template waveform, while faithfulness is the overlap maximised over only 
the extrinsic parameters of the template. Faithfulness is a measure of how good the template 
waveform is in both detecting a signal and estimating its parameters. However, the fitting 
factor is aimed at finding whether or not a template family is good enough in detecting a 
signal without reference to its use in estimating the parameters. 

The fitting factors and faithfulness of three different template families, using the hybrid 
waveforms as target signals, are plotted in Figure [2j Maximization over the intrinsic 
parameters is performed with the aid of the Nelder-Mead downhill simplex algorithm. Low 
frequency cutoff is chosen to be equal to 40Hz for Initial LIGO and 20Hz for Virgo and 
Advanced LIGO. As expected, PN inspiral templates produce very good overlaps with the 
target signals in the low-mass regime (where inspiral is the dominant part), and ring-down 
templates produce good overlaps in the high-mass regime (where ring down is the dominant 
part). The black-hole coalescence templates continue to produce very good overlaps over 
the entire mass range. Note that, if we assume homogeneous and isotropic distribution of 
sources, the fraction of sources detectable by a template family is proportional to the cube 
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of the fitting factor [39]. The 'event loss' due to mismatch between the template and the 
true signal is larger than the canonical 10% when the fitting factor is less than 0.965 fjj] 
This figure suggests that the template family proposed in Section [3] can be used to search 
for binary black-hole coalescences over almost the entire mass range considered here losing 
no more than 10% of the events that are detectable by optimal filtering. Binaries with 
total mass < 15M Q are detectable using PN inspiral templates with < 10% event loss, 
while the binaries in the mass range considered here cannot be detected with < 10% event 
loss using black-hole ring-down templates |JJ Faithfulness of the coalescence templates is 
also almost always greater than 0.965 (comparable to the fitting factors), while that of the 
other templates are considerably smaller in general. This also means that the parameters 
estimated by the PN and ring-down templates will be biased significantly. This will be 
studied in detail in a forthcoming work [33, ■ 

We can also calculate the 'distance reach' of these searches. Since the template 
waveforms described in Eqs.([5]) and ^ are shown to be very close (fitting factors > 0.95) 
to the hybrid waveforms, the effective distance c? op t to binaries producing a certain optimal 
signal-to-noise ratio (SNR) p at the detector can be computed analytically using these 
template waveforms, as: 

^c ff (/)d/] 1/2 _ (8) 



d Q pt 



s h (f) 



Since the fitting factor (FF) is the fraction of the optimal SNR that can be achieved 
using a sub-optimal filter, the effective distance d su bopt to the (optimally-oriented) binaries 
producing a sub-optimal SNR p by a template family is given by d su bopt = d opt FF. Figure[3] 
compares the effective distance to optimally-located and- oriented binaries producing a sub- 
optimal SNR of 8 at the detector output using the three different template families discussed 
above. For the black-hole coalescence templates, the horizon distance reaches peak values 
of around 760 Mpc (at 150 M ), 950 Mpc (325 M ) and 13.3 Gpc (225 M Q ) for Initial 
LIGO, Virgo and Advanced LIGO, respectively. For PN inspiral templates, the peak values 
are 630 Mpc (160 M Q ), 770 Mpc (325 M ) and 9.2 Gpc (250 M ), while for the ring-down 
templates, the corresponding values are 640 Mpc (150 M ), 780 Mpc (325 M ) and 10.6 
Gpc (225 Mq). It may be noted that, since the event rate is proportional to the cube of 
the distance reach, a 20% loss in the distance reach means a 50% loss in the event rate. 



5. Summary and future work 



Recent progress in the theoretical modelling of coalescing binary black holes has important 
applications in the search for gravitational waves from binary black-hole coalescences. 
'Complete' gravitational waveforms can be constructed by combining results from analytical 
and numerical calculations. These waveforms can be parametrised to produce analytical 
waveform templates which can be used to densely cover the parameter space of the binary 
that will be searched over by matched- filtering techniques. This template family will allow us 
to coherently search over all the three stages of the binary. The advantages of this 'coherent 
search' include improved SNR, and hence improved distance reach for the search, potential 

|| It should be noted that fitting factor is not the only consideration in an actual search strategy. Other 
factors such as computational cost and false alarm rate also play a decisive role in the choice of a template 
family. However such a detailed study is out of the scope of this paper. 

1 Note that the overlaps are maximised over the initial phase ipg of the ring down also, unlike what is 
proposed in Ref. |36| . 
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Figure 3. Left plots shows the horizon distances of PN inspiral, ring down and 
coalescence templates in the case of Initial LIGO (thin lines) and Virgo (thick lines) 
noise spectra. Right plots show the same for Advanced LIGO. Dots correspond to 
PN inspiral templates, triangles to ring-down templates, and squares to coalescence 
templates. Horizontal axes report the total mass of the binary, and the vertical axes 
report the effective distance to optimally-oriented equal-mass binaries producing a sub- 
optimal SNR of 8 at the detector output. The sharp drop in the PN horizon distance is 
a result of the (different) lower cutoff frequencies of the detectors. 



reduction of the false-alarm rate and improved parameter estimation. The application of 
the complete waveforms is not limited to template-based searches. For example, these 
waveforms can also be used in the burst searches, customised for detection of signals from 
coalescing black-hole binaries. In this case, a small bank of representative waveforms can 
be used to survey the parameter space of the binary |40j . The complete waveforms can also 
be used to estimate the efficiency of different search methods. A preliminary comparison 
of three different template-based searches is presented in this paper. More robust ways of 
comparing the efficiency of different searches is an ongoing effort |34j . 
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